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Abstract 

Recently there has been a growing interest in nanodroplet engi- 
neering and the dynamics of vaporization and condensation of fluids 
at the nanoscale. The applications of nanodroplet engineering is in 
the biomolecular field, physical chemistry, and engineering power pro- 
duction and storage devices technology... other applications are in 
aerosols, fuels and turbulence in power source design. In this paper a 
theoretical model for multiphase nanofluids is derived from informa- 
tion theory and temperature dependent distributions are obtained for 
the Lennard- Jones potential interactions of the nanofluid droplets. 

The simple paradigm of a nanofluid is an inert material such as Argon, as 
a fluid at low temperatures, in a droplet of a diameter a few nano meters in 
scale. Such a nanodroplet contains from the tens to the hundreds of Argon 
atoms in a fluid phase. The kinematics of the fluid nanodroplet as it is 
accelerated and it interacts with other materials, gas, liquid, molecules or 
solid state materials is the dynamics we wish to address here. 

A simple example for what we will model by theory here is the nan- 
odroplet of Argon, say containing a few hundred Argon atoms, moving at a 
uniform velocity towards a non chemically reactive solid surface that could 
represent in more complex experiments or applications a reactive or non re- 
active surface of a device for power production such as in fuel cell walls, 
a porous membrane in biological cellular reactions, or reactive chemistry 
analysis such as lab-on-a-chip. The medium here is a vacuum though more 
complex arrangements could be of secondary gases or liquids. 

At the nanoscale the aggregation of atoms and the screening effect that 
occurs as the number of atoms increases towards a size of the micron range 
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is less pronounced , and quantum mechanical effects tend to emerge and 
be more pronounced. This is a general phenomenon and is noted in large 
molecules as opposed to assemblages of molecular chains, in nanoparticles as 
opposed to bulk solids. Examples are the magnetic hysteresis in Gold and 
Palladium nanoparticles that are diamagnetic in Bulk [TJ . Another bulk 
equivalently large scale effect are the continuum bands of electronic transport 
and the breakdown of these bands and with the subsequent localization of 
electronic orbitals described by hopping Hamiltonian methods. 

An observed effect at the nanoscale of the breakdown of continuum elec- 
tronic bands is the quenching of plasmon oscillations. For a discussion of 
these quantum nano scale effects see any recent book on nanoparticle engi- 
neering. Electronic Plasma is a Fermi fluid description of the collective bulk 
effect of interacting electrons in continuum bands. This has been used to very 
good design advantage in absorption and emission spectrum of interaction 
with photons and the electronic plasma in image enhancement and medical 
applications. 

The emergent quantum phenomenon in nanofluid droplets is expected to 
be similarly novel and of interest in engineering and applications. There is 
then great interest in accurate models of multiphase processes and the dy- 
namics of nanofluids in general and of nanodroplets in particular. Multiphase 
processes of phase transitions are of importance as discussed in vaporization 
applications such as in modeling vaporization in turbulent regimes and fuel 
mixing, in biomolecular applications, etc. 

Theoretically the modelling of nanoscale fluids should be done at the fully 
quantum mechanical level. This means the interacting many-particle physics 
approach, many particle Hamiltonians and diagrammatic methods. However 
the screening effect of aggregate atoms allows for approximate descriptions 
at various length nano and mesoscopic scales where either fully quantum 
mechanical effects emerge or partially screened effects are noted and can be 
described by semi classical approximations at the longer length scales from 
the nanometer 10~ 9 m to 10 _6 m micrometer range. At larger length scales 
it is expected that classical descriptions of fluids are accurate for 'bulk' fluid 
description. 

1 Fully Quantum Mechanical Theory 

A fully quantum mechanical theory of nanofluid droplets has already been 
written and is published in the many particle theory of Fermi systems phase 
transitions. The review articles by Mattuck and Johanssen are of very in- 
fluential importance in this theoretical approach [21 EJJ. The theory is a 



2 



description of a fluid and a vapor phase coexisting and described simultane- 
ously in a matrix propagator or Green's function. The quantum mechanical 
Hamiltonian is a many-particle Hamiltonian of a form that can be written in 
occupation number formalism or in second quantization as sums of creation 
c' and annihilation c operators of the observables of kinetic and potential 
energies. 

H = J2 € k C h C k+ V k,l,m,n cA[ k C f Wn- (1) 

k k,l,m,n 

The multiphases of interest are solid, liquid and vapor or gaseous phases. 
These are described by the inclusion of translation periodic vectors for the 
solid phase such that the Hamiltonian of Eq.([T]) becomes periodic with re- 
spect to the lattice structure of the solid phase, the lattice periodicity vector 
we set to be a constant Q Q , Q' Q for simplicity. The Hamiltonian of the solid 
phase becomes 

inf inf 

H = S ^ t K+Qo C ^k+Qo C k + ^k+Q X,m+Q^ ,n C ^k+Q C l C ^rn+Q' C n-- ( 2 ) 

k,Q k,l,m,n,Q ,Q' 

As the temperature is increased, the solid melts to a fluid phase, and the 
Hamiltonian is described by a similar Hamiltonian, however the lattice sym- 
metry or periodicity is broken, and the lattice vector which was set to be 
a constant Q a for simplicity as in Eq.(j2]) is now a nonperiodic momentum 
vector and the fluid is described by a Hamiltonian of the form of Eq.([T]). 
Note that the solid is described by an effective mass electronic Hamiltonian 
describing band electronic effective mass. This effective mass simple descrip- 
tion is no longer applicable in the fluid phase and actual electronic mass or 
effective atomic charge mass and effective nuclei is to be used in the Hamilto- 
nian approximation, or a full interacting Hamiltonian of the protons and the 
electrons with description of the interactions of the protons and the electrons 
in each atom of the fluid and between atoms.... This would mean adding in 
superposition to Eq.(TfJ) terms of the potential interaction between protons 
and terms of potential interaction between protons and electrons such that 

< KE P > + < V p - P >= ^eVW+ Yl yP ^ '^rC 1 W 1 <fr (3) 

o o,p,q,r 

and interaction terms between the protons and electrons of the similar form 
< Vp-e >= Yl V p ~ e 3,p,q,rC* oCpf^ qCr, though the choice of potential can differ 

o,p,q,f 
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for effective charge interactions in band material, for effective Hydrogen like 
atoms, etc. 

The gaseous or vapor phase is again described by a similar Hamiltonian 
as in Eq. (TTJ , the difference here being that again as in the fluid phase Hamil- 
tonian we can no longer use the approximations we used in the solid phase 
of effective mass band electrons. Also, the Hydrogen like atom with effective 
charge and effective nuclei that is useful for simple atoms description in the 
fluid phase is applicable for simple atoms. Otherwise we would again have 
to fully describe each electron and proton in each atom in the fully interact- 
ing picture we have discussed, with electron and proton and proton-electron 
interaction potentials in superposition with the electronic Hamiltonian. The 
finite temperature and the number of the atoms in each phase are the impor- 
tant considerations we discuss next. The number of atoms in each phase is 
dictated by the chemical potential; the number of particles is the summation 
of the Fermi distribution function as in < iV >= e p(* k -n) +1 ■ How one par- 

k 

titions the groups of electrons , protons, and atoms they comprise depends 
on the level of approximation used as discussed. 

The Hamiltonians can be written as a matrix Hamiltonian, or as one 
single Hamiltonian using spinors. The matrix Hamiltonian also means that 
the Green's functions of the Hamiltonians can be written as a matrix Green's 
function. The Hamiltonians would then be written as 



< H >-- 



^ H solid— Solid H solid— Liquid H Solid— Gas / 'vapor 

^ H Liquid— S olid H Liquid— Liquid H Liquid— Gas/ 'vapor ■ (^0 

^-Gas/vapor— Solid ^ Gas j vapor— Liquid T Gas / 'vapor— Gas / 'vapor 



The Green's function in quantum mechanics is simply the response func- 
tion, and it can be described as the inverse of the Hamiltonian as G r (k, oS) = 
t — here the superscript G r indicates the retarded Green's function and 
the imaginary infinitesimal +ir] in the denominator insures no divergence in 
the Green's function by shifting the argument into the complex plane and 
the total kinetic and self energy is written in shorthand as v? = — iH{k, u) 
[I] . Physically the causality of incoming or outgoing waves is represented by 
this shift into the complex plane. 

The Green's function for the total Hamiltonian of the coexisting phases 
is then also a matrix, and can be written as in Eq.fllJ), 



G 



Qs—s 9s— L 9s— v 

9l-s 9l-l 9l-v 

9v—s 9v—L 9v—v 



(5) 



and here the Green's functions are subscripted in simple notation such as the 
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solid-vapor interacting Green's function g s _ v and so on for the solid-solid , 
liquid-liquid and the vapor-vapor (gaseous) phases. 

The observables of interest are the total number of particles, the number 
of particles or atoms in a particular phase, the rate of vaporization or of con- 
densation or solidification, temperature and energy dependences of the phase 
transitions. These are obtained in expectation values of the Green's func- 
tions, also to be considered as probability functions analogues here, in sum- 
mation with the observable of interest. The total number of Argon atoms for 
example, < iV >=< N s - a > + < Nl-l > + < N v - V > is obtained by the di- 
agonal trace of the Green's function matrix since < iV >= J G(k, lo) ,£ 3 ^i 
and the trace consists of the individual phases such as for example < N s _ s >= 
J g s - s (k, lo) Fourier transforming from the energy to the time domain 

and taking a time derivative can describe vaporization from the liquid phase 
for example 9<N Q t ~ v> = ~§i f 9L-v(k, t) o^ks dt, while a similar transform and 
time derivative can describe sublimation from the solid to the vapor phase 

a<N s -v> _ a r tt f \ dk j f 
at ~ at J ys-v^, ^ (211)3 

The solution of the equations for the observables, as they are rather com- 
plicated mathematical expressions, can be done by computer calculations 
of approximate summation methods such as the Hartree-Fock approxima- 
tion, mean field and self consistent methods and dynamic mean field and 
other similar summation methods, and only rarely in the formal diagram- 
matic theory are closed form solutions possible, as in linearized mean field 
like approximations, other approaches to analytic solutions are possible by 
transformation methods, and these are discussed in later work. 



2 Semi Classical Theory. Lennard Jones Flu- 
ids. 

The nanofluid droplets have been described reasonably successfully in the 
literature by semi classical methods [5j UJ. The approaches have been to 
treat the nanofluid droplet by hydrodynamic means, the vapor by Boltzmann 
equation, and to evolve the number of particles and obtain rates of change in 
number of particles. Approximations are made such as Boltzmann equation 
[7J collision terms of a certain order, an example is the force collision terms 
up to a third order expansion by [6]. 

In hydrodynamics approaches, full analytic solutions are not obtained to 
the hydrodynamic equations... approximations of linearization and similar 
assumptions for solution are made such as in [5] , examples of approximations 
are of viscosity and compressibility. 
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Another approach is possible. Given the success of recent molecular dy- 
namics simulations in describing experimental results, approximating the 
nanofluid as a Lennard- Jones fluid seems reasonable. Then a similar the- 
ory can be made at the Lennard- Jones approximation level that is applicable 
in the range of applicability of the Boltzmann and hydrodynamics equations' 
scale of applicability. That is from the hundreds of atoms to the micrometer 
scale. From an analogy to the fully quantum mechanical phase transition 
theory discussed in the previous section, a theory of coexisting phases of 
nanofluid droplets can be obtained. The advantages would be the simulta- 
neous description of co-existing phases by analogy to the matrix formulation 
discussed, long range order parameter approach as in the diagonal trace and 
summation to obtain particle or number of atoms in the individual phases. 

A theoretical advantage of the full quantum theory is the Hamiltonian su- 
perposition that allows us to include external perturbations and interaction 
terms. From the point of view of the approach we will pursue in order to de- 
rive a Lennard- Jones fluids coexisting phases theory, the superposition prin- 
ciple is replaced by the Legendre transform of the maximum entropy method 
of theory derivation, and equivalently the information theoretic method. Al- 
ternatively a Lennard- Jones semi classical Lagrangian or Hamiltonian may 
be used, and dynamics with or without a heat bath and therefore Fokker- 
Planck equations obtained in matrix form. We proceed from the information 
theoretic and equivalently maximum entropy derived theory approach. 

2.1 Derivation 

The entropy measure to be used we choose to be the Gibbs-Boltzmann en- 
tropy, < S >= — c J P(r, t)lnP(f, t)df which can be generalized to the nonex- 
tensive entropy of C. Tsallis [TU|ITTj for better accuracy and other novel effects 
such as the q-parametrization of the deviation of the PDF from the Gaussian 
to the power-law PDF. Note that nonlinearity due to interactions naturally 
or physically deviate the PDF from the Gaussian distribution towards the 
power-law distribution. The nonextensive statistics has been noted for theo- 
retically producing this form of power-law statistics from a formally similar 
mathematics to the traditional thermodynamics and information theory. The 
connections are more profound and review articles and applications can be 
found at [10J . 

The maximum entropy approach is then of the maximization of the en- 
tropy measure constrained with observables of the system in a Legendre 
transform with Lagrange multipliers setting the weighted units to the en- 
tropic measure's. 

The choice of observables is one of choosing the minimum set of observ- 
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ables that describe the physics of interest. These observables are described 
mathematically as functions of the other observables chosen as representative 
of the observable dynamics of the system. A discussion of this was made by 
the original author of maximum entropy theory [HJ [9] E. T. Jaynes. 

We choose observables analogous to the quantum mechanical observables 
in the fully quantum mechanical theory. The momentum or kinetic energy ob- 
servables of quantum mechanics correspond to the first and second moments 
of the variables of position, the potential energy or interaction potential 
terms are similarly included here though we choose the Lennard- Jones po- 
tential and not the Coulomb potential, and the addition of velocity variables 
insures that a connection to classical equations of motion such as hydrody- 
namic and Boltzmann is made. The observable moments are: 

N v 

M 2 (Axij) = ^2< (A%- < Axij >) 2 > 
= J (Axij— < Ax^ >) 2 P(Ax~ij, Av~ij, At)dAxijdAv~ij (6) 

N v 

M 2 {Av lj ) = < < Avij >f > 

hi 

= j(M„- < A % »V(A % Av^, At)dAxijdAvij (7) 



N v 

^ v vapor -vapor ^ ^ \ g 12 ) 

N v 

= E / - ^) P ( A % A %-> Af)dAa$dAu£ (8) 

i,j ' '"' J '''I 

Here we have chosen to use the Lennard- Jones potential in full form as ap- 
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propriate to the level of mechanics and scale. Note that in solid form a 
harmonic oscillator approximation has been used by some researchers and in 
fluid phase this is somewhat questionable. The liquid and the vapor phases 
interactions are chosen to be described by Lennard- Jones potentials. The 
Lennard- Jones potential is a scale dependent approximation for screened in- 
teractions of coulomb like potential and include the repulsive term where no 
two atoms can occupy the same space. The potential behaves as an inverse 
Coulomb-like potential at far enough distances. In quantum mechanics 
the repulsion of atoms at almost vanishing distances is the Fermi exclusion 
principle statistics or for zero or integer spin atoms and particles the Bose 
statistics which as known can overlap and condense. As discussed in the pre- 
vious sections the potential in the quantum mechanical interaction theory is 
a coulomb or screened coulomb potential. 
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The maximization of the entropy measure with the observables constrain- 
ing the maximization allows us to derive a least biased distribution [8J. As 
an example we derive the vapor phase least biased distribution, 



5 < S > +5{(3M 2 (Ax ij ) + aM 2 {Axij) + 7 < V v - V >] = 0. 



(9) 



The least biased distribution is 



/3(Ax ii -<Ax ij >) 2 



■a(At/ i:) -<AC ii >) 2 -7( 




P{Axij,Av ih At) 



e 



(10) 



Z(At) 



The Lagrange multipliers can be evaluated from the identity 



1 dZjt) 
Z(t) dp 



< (Ax^- < Axij >) 2 >, 



(11) 



similar relations between observables and for the other Lagrange multipliers 
can be evaluated. 

The partition function is the inverse of the normalization of the proba- 
bility distribution function and as usual is the sum over the distribution and 
its variables as in Z(t) = f P(x,y, ..,v z ,t)dxdy..dv z . 

The matrix of co existing phases here is obtained similarly to the fully 
quantum mechanical theory. The method of matrix or spinor Hamiltonian 
and derivation of Green's functions is discussed in the quantum mechanical 
case and will not be revisited here, we proceed to the discussion of the vapor 
and liquid phases and the interaction terms. 

Due to the use of similar Lennard- Jones interaction potentials between 
the liquid-liquid atoms, the liquid-vapor atoms, and the vapor-vapor atoms 
and the solid-solid atoms, the mathematical distinction between the vapor 
and liquid phases is one of number of atoms in an individual phase and 
not the analytic form of the least biased distribution function of particles 
in a liquid solid or vapor phase. This is made clear when we note that the 
nanoffuid droplet of number of atoms < Nl >, where < N >=< Nl > 
+ < N v > + < N s > is a constant number and is obtained by summing the 
distribution of atoms over the variables of the observables 
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< N L >= ^2 / P(&Xij, A%, At)dAxijdAvij (12) 

N v . 

<N V >=^2 / P(Ax i:j , Avij, AtfdAxijdAvij, (13) 



<jV M >= PiAxij^Vi^A^dAxijdAvij, (14) 



; (is) 

The use of the Lennard- Jones potential in both the vapor and fluid phases 
is the cause of the formal similarity of the terms here. If we had chosen to use 
the harmonic oscillator potential for the solid phase interactions, for one rea- 
son or another, and the Lennard- Jones potential for the vapor and a similar 
screened Coulomb potential for the liquid phase, or had chosen to treat the 
solid phase as fully quantum mechanical while the liquid phase was treated 
as a Lennard- Jones liquid while say the vapor phase was treated also by a 
Lennard- Jones or by a Boltzmann equation PDF, the mathematical analytic 
similarity would no longer exist, and the off diagonal terms in the matrix 
of the Argon atoms distribution functions would no longer seem redundant. 
The matrix of these distribution functions is in direct analogy to the fully 
quantum mechanical theory 



Ps—s Ps—L Ps—v 

p L -s p L - L p L - v 

Pv—s Pv—L Pv—v 



(16) 



Other instances where offdiagonal terms and therefore offdiagonal long 
range order gives no information is the quantum mechanical mathematical 
analogue of the mean field Ferromagnetic model, where the preferred di- 
rection and the linearization of the field terms makes the offdiagonal terms 
redundant. As such they are constant terms in the original Hamiltonian 
superposition and are zero terms in the Green's function matrix analogous 
to the PDF matrix here, which then contains only diagonal terms and only 
diagonal long range order. 
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Here then the analogy to the diagonal long range order distributions (re: 
Green's functions) in the quantum mechanical phase transition theory has 
the diagonal PDFs as the distributions importance in determining the observ- 
ables of interest such as the total particle number and inter-relationship of 
the phases in determining rates of change of one phase's number of particles 
or atoms as dependent upon the other phases' number of atoms or particles. 
The redundant seeming offdiagonal terms are also similarly able to provide 
the information about the rates of change of number of atoms from one phase 
to another, and this is discussed subsequently. Also to be considered is the 
symmetry between offdiagonal terms such as p s -L <=? => Pl-s- 

The similarity of potential terms used in all the phases, here the Lennard- 
Jones potential used in all three co-existing phases makes the coexisting va- 
por and fluid and solid terms of individual phases described by the diagonal 
terms, the analogues of the diagonal long range order in the quantum me- 
chanical phase transition theory, similar in form to the offdiagonal terms and 
denotes the importance of the number of atoms or particles in each phase 
as the distinguishing feature in our simplified approach. We note again that 
the use of dissimilar potentials or even statistics such as quantum mechani- 
cal descriptions of one phase and semi-classical descriptions of another phase, 
would result in distinct differences in analytical form of the distributions and 
the perceived redundancy of offdiagonal terms would no longer exist, and it 
is possible that the importance then of offdiagonal long range order would 
outweigh the diagonal terms. Note that we use the traditional terminology 
of diagonal and offdiagonal long range order in analogy to the fully quantum 
mechanical phase transition theory. 

The similarity of potential terms used in all the phases, here the Lennard- 
Jones potential used in all three co-existing phases makes the coexisting va- 
por and fluid and solid terms of individual phases described by the diagonal 
terms, the analogues of the diagonal long range order in the quantum me- 
chanical phase transition theory, similar in form to the offdiagonal terms and 
denotes the importance of the number of atoms or particles in each phase 
as the distinguishing feature in our simplified approach. We note again that 
the use of dissimilar potentials or even statistics such as quantum mechani- 
cal descriptions of one phase and semi-classical descriptions of another phase, 
would result in distinct differences in analytical form of the distributions and 
the perceived redundancy of offdiagonal terms would no longer exist, and it 
is possible that the importance then of offdiagonal long range order would 
outweigh the diagonal terms. Note that we use the traditional terminology 
of diagonal and offdiagonal long range order in analogy to the fully quantum 
mechanical phase transition theory. 

The total number of particles is a constant in our simple 100 atom Argon 
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nanodroplet as pointed out, and is the sum described as < N >=< N L > 
+ < N v > + < N S >. The number of liquid particles at a time t is the 
< Nl >=< N > — < N v > — < N s >, and the rate of change of this number 
is the partial time derivative d< Q t L> = ~ 9< Q t v> ~ 9< dt > m This corresponds to 
the time rate of change in the number of liquid atoms that are being vaporized 
or solidified. Other interaction distributions, say between the liquid and 

N V ,N L 

vapor atoms < N v -l >= J P(Axij, Avij, At)dAxijdAvij are useful in 

ij 

rates of transition, where a time derivative is the rate at which liquid Argon 
droplet atoms vaporize to the gaseous phase. These offdiagonal terms are 
important in offdiagonal long range order derivation, and here in rates of 
transition from one phase to another. They are the semi classical analogues 
of the quantum mechanical Green's functions of interacting vapor-liquid and 
liquid-vapor and solid-liquid and solid-vapor Hamiltonians. Differences are 
that the quantum mechanical off diagonal Green's functions are Hermitian 
transposes of each other. Here the analogue offdiagonal PDFs are to be 
considered from the point of view of reversibility and mathematical similarity. 
Transpose of coordinates in our two-point PDF functions does not result in 
a change of sign for de-meaned second moments, however it is possible that 
with moments about the mean, the mean or drift terms could introduce 
a sign change upon interchange of coordinates. Care must be take in any 
assumptions about the symmetry of offdiagonal terms. The even power form 
of the Lennard- Jones potential implies that interchanging the interaction 
coordinates does not cause a change in sign overall however this might not be 
the case for dissimilar potentials and statistics between phases as mentioned. 



2.2 Derivation of Solutions 

An interesting point to note is that we can obtain closed form solutions for 
the < N L > and < N v > and < N s > PDFs . The PDF for the nanofluid 
droplet is 



-/3(A^ 3 -<A^ 3 » 2 - a (A^-<A^» 2 - 7 (^-2 ¥f ) 

P(Ax tJ , Av^, At) = - _ (17) 

The Lagrange multipliers are evaluated as discussed from the relation- 
ship between partition function and the observables = — M 2 (Axij) and 
similarly for the other observables. Note that the Lagrange multipliers of the 
second moments are related to the diffusion coefficients by the temperature- 
diffusion coefficient thermodynamics relation D — js, and < (:%) 2 ) >= 2Dt. 
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The analytic solution of interest here is the one of the transformation of the 
Lennard- Jones potential to a drift coefficient term. This is done by uncom- 
pleting the square of the Gaussian 

p(X) = exp(X 2 -2X <X> + <X > 2 )/2Dt (18) 

and comparing it with the dependent terms of the de-mean (or constant 
or linear and renormalized mean) Gaussian 

P(r iv ....) = e- /3 «4- 2 ^<^>+<^> 2 )-^ J (^)) e e (19) 

The Lennard- Jones potential terms can be seen to be terms that can be 
transformed to a drift coefficient or the first moment by un-completing the 
square. We obtain 



< Tij >= r i:j + Jrij 2 + -LJ( rij ). (20) 



2 

r 

2.3 Analytic distribution functions 

With this the dependent terms become a simple e~Pi r io-< r n>)' 2 Gaus- 
sian which by the identity = —M 2 (Axij) gives the simple diffusion- 
temperature relation = 7^-7 with the temperature from thermodynamics 
= — 

The resulting PDF is a Gaussian with the coordinates expressed in spher- 
ical coordinates for simplicity. The integration for the partition function or 
normalization can be accomplished, however we use the relationship (3Z(t) 2 = 
2n and the resulting distribution is 

-(^,-<^,» 2 -^-^li^ 2 

g 2D r t g 2D v rt 

P( rij , ....) = 3 3-. (21) 

y/4:TlD r t yjAllD v rt 

Note that aside from the terms which have the complicated Lennard- Jones 
potential dependence and its radial coordinates dependence Eq. (120]) added 
to the radial coordinate linear drift, the other variables such as the velocity 
magnitudes in different coordinates are all simply zero mean < 0^ >= 0, 
constant drift < 6ij >= b or linear drift terms corresponding to accelerated 
motion < 6{j >= c6ij + 6; Note that a constant drift value mean is useful 
in uniform motion descriptions, and a linear < 6ij >= c8ij + b drift term 
corresponds to a constant driving force term. 
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3 Computation and numerical simulation. 



The simulation of the statistical dynamics of a nanofluid droplets is accom- 
plished by using these probability distribution functions in a Monte Carlo 
or equivalently Stochastic simulation. In the next of the series of papers on 
nanofluid droplets, the simulation of nanofluid Argon droplet is performed 
using a particular set of parameters. We choose to simulate a nano particle 
of number of Argon atoms of 100, at room temperature 300K, in a vacuum 
and with no surfaces to interact with for the initial simulation of vaporiza- 
tion. The simple high temperature droplet parameters mean we can neglect 
the solid phase and the 3X3 matrices become 2X2 matrices with no loss of 
generality. The Diffusion coefficients are obtained from the inverse of the 
temperature and the Boltzmann constant as discussed. The Lennard- Jones 
parameters are obtained from CHARMM Molecular Dynamics Argon file 
data from experimental measured potential parameters. The trajectory of 
< N >= 100 Ar atoms is obtained from the stochastic differential equations 
of the individual atoms. The initial fluid configuration is assumed to be a 
fixed configuration, and initially the fluid droplet is at zero coordinate. The 
initial configuration is assumed to be almost solid Argon which is FCC face 
center cubic, and the average nearest neighbour distance is 3.6 Angstroms, 
and the next nearest neighbour distance is obtained from that by geometrical 
analysis and so on. Alternatively we simulate a periodic initial configuration 
with variation of nearets neighbor distances about the mean nearest neighbor 
distance of the FCC solid phase value. This initial configuration is almost 
immediately melted as the high room temperature in comparison to fluid 
Argon temperature of approximately 67K causes the atoms to move apart 
rapidly. The movement is somewhat restricted by the potential which is at- 
tractive at long distances and repulsive at short distances, and the fluid phase 
begins to vaporize as atoms separate from the nano fluid droplet and move 
independently off. As in other literature, we define the vaporized atoms as 
atoms that move 3 nearest neighbor distances away from the nearest fluid 
atoms that are all within 1.5 times the nearest neighbor distances to each 
other. In simulation coordinates, we can plot the trajectories of atoms and 
can observe the vaporization of Argon atoms from the nanofluid droplet. The 
simulation is run until half of the atoms are vaporized. We also estimate this 
by tracking the rate of change of number of fluid atoms. 
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4 Conclusion 



In this article we have derived a Lennard- Jones fluid theory of nanofluid 
droplet statistical dynamics from an information theory and equivalently 
maximum entropy approach. We have derived analogously to the fully quan- 
tum mechanical theory of phase transitions of fluids a matrix of probability 
distributions and therefore number of Argon atoms that are in a particu- 
lar solid, fluid or vapor phase at any given time. We discuss the current 
approach's analogues of diagonal and off diagonal long range order in par- 
ticle or atom number summation and rates of transition from one phase to 
another. We then solve analytically a particular distribution function and 
obtain analytic forms that are normalized. This analytic form is due to the 
radial dependence of the Lennard- Jones potential. We then discuss research 
of computational simulations performed on nanofluid droplets which we will 
report in detail on for materials of Argon for simplicity, water as a check of 
more complex nanofluids and for future applications to biomolecular model- 
ing in nanofluids. 
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